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1. INTRODUCTION 


1.1 Rotor Blade Design 

The design of helicopter rotor blades Involves not only considerations of 
strength! survivability, fatigue, and cost, but also requires that blade 
natural frequencies be significantly separated from the fundamental aerodynamic 
forcing frequencies (e.g. Ref. 1). A proper placement of blade frequencies 
is a difficult task for several reasons. First, there are many forcing 
frequencies (at all integer-multiples of the rotor RPM) which occur at rather 
closely-spaced intervals. For example, 5/rev and 6/rev are less than 20% 
apart. Second, the rotor RPM may vary over a significant range throught the 
flight envelope, thus reducing even further the area of acceptable natural 
frequencies. Third, the natural modes of the rotor blade are often coupled 
because of pitch angle, blade twist, offset between the mass center and 
elastic axis, and large aerodynamic damping. These couplings complicate the 
calculation of natural frequencies. In fact, the dependence on pitch angle 
makes frequencies a function of loading condition, since loading affects 
collective pitch. Fourth, the centrifugal stiffness often dominates the 
lower modes, making it difficult to alter frequencies by simple changes in 
stiffness. 

In the early stages of the development of the helicopter, it was believed 
that helicopter vibrations could be reduced (and even eliminated) by the correct 
choice of structural coupling and mass stiffness distributions. However, it 
is easy to imagine how difficult it is to find just the proper parameters such 
that the desired natural frequencies can be obtained. The difficulties in 
placement of natural frequencies have led, in many cases, to preliminary 
designs which ignore frequency placement. Then, after the structure is 
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’’finalized" (cither on paper or in a prototype blade) » the frequenciee are 
calculated (or meaeured) and final adjuatmenta nadct Reference (2) deacribea 
the development of. the XH-17 helipopter in which a 300-lb weight waa added to 
each blade in order to change the spanwiae and chordwiae maaa diatribution and 
thereby move the firat flapwiae frequency away from 3/rev. However* theae 
typea of alterationa are detrimental to blade weight* aircraft development 
time* and blade coat. In addition* correctiona uaually are not aatia- 
factory, and the helicopter is often left with a noticeable vibration problem. 

The state-of-the-art in helicopter technology is now to the point* 
however* that it should be possible to correctly place rotor frequencies 
during preliminary design stages. There are several reasons for this. First* 
helicopter rotor blades for both main rotors and tall rotors are now being 
fabricated from composite materials (Refs. 3 and 4). This implies that the 
designer can choose* with limited restrictions* the exact El distribution 
desired. Furthermore* the lightness of composite blades for the main rotor 
usually necessitates the addition of weight to give sufficient autorotatlonal 
blade Inertia. Thus* there is a considerable amount of flexibility as to how 
this weight may be distributed, gecond* the methods of structural optimization 
and parameter identification are now refined to the point where they can be 
efficiently applied to the blade structure. Some elementary techniques have 
already been used for the design of rotor fuselages (Ref. 5). It follows that 
the time is right for the use of structural optimization in helicopter blade 
design. Some work on this is already under development, and* although not 
published, some companies are already experimenting with the optimum way to 
add weight to an existing blade in order to Improve vibrations. 
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Th« purpo 9 « of th« work diacuoaod h«r« la to invaatlgi a tha poaalbilltlaa 
<aa wall aa tha llmi tat Iona) of tailoring blada naaa and aciffnaaa diatrlbutiona 
to g:;va an optlnuan blada daaign in tarma of waight, inartia, and dynamic 
characteriatica. 

The major objactivea of tha work are: 

1) To datannina to what extant changea in maaa or atiffnaaa diatributlon 
can be uaad to place rotor frequanclea at deaired locatlona. 

2) To eatabllah theorfitlcal limits to the amount of frequency shift. 

3) To formulate realistic constraints on blade properties based on weight » 
mass moment of inertia, size, strength, and stability* 

4) To determine to what extent the hub loads can be minimized by proper 
choice of EX distribution* 

5) To determine if the design for minimum hub loads can be approximated 
by a design for a given set of natural frequencies. 

6) To determine to what extent aerodynamic couplings might affect the 
optimum blade design. 

7) To determine the relative effectiveness of mass and stiffness distribution 
on the optimization procedure. 

8) To determine to what extent an existing blade could be optimized with 
minimal changes in blade structure. 

9) To develop several "optimum profiles" for rotor blades operating under 
various standard conditions. 

The work is. to focus on configurations that are simple enough to yiuld 
clear, fundamental insights into the structural mechanisms but which are 
sufficiently complex to result in a rt^allstic result for an optimum rotor blade. 
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X.2 Qvarvittw o£ Optlioal Structural Daatgn 

Moat approachaa to optimal atructural daaign may ba claaai£lad into thraa 
catagoriea. (For recent review articlaa aee Ra£a» 6 and 7.) One avxh category 
is *Srariational methoda." Theae generally rely on techniquea from the mathe- 
matical theory of the calculua of variationa, and, when applicable, often 
provide uaeful phyalcal insight into the nature of an optimal design • 
Unfortunately, only relatively simple problems can be solved by this approach, 
since the matltematics becomes intractable when complex engineering structures 
are considered* 

A second category of atructural optimization techniques consists of the 
application of mathematical programming methods together with the discreti- 
zation of the structure by finite element techniques. This approach to 
optimization was founded in I960 (Ref. 8) with the hope that more complex 
structures could be analyzed than were possible when using the analytical 
techniques of the calculus of variations. However, in the late 60* s It 
became apparent that mathematical programming methods had limitations of their 
own, namely, unacceptably long computation times occurring when the number of 
design variables become large (over 20-100, depending on the type of structure). 
Fortunately, several Improvements developed over the last few years appear to 
have significantly extended the capability of the mathematical programming 
approach, and, as a result, it is this approach we intend to draw upon for 
solution techniques in the proposed research. In a later section of this 
proposal, after our design problem has been formulated, we will discuss these 
recent improvements in the approach. 

A third cn'-egory of structural optimization approaches is the ’’optimality 
criterion” approach in which an equation expressing some necessary condition of 
optimality is used as the basis for constructing an iterative (successive 


rfi-dttsign) procsdur*. OrlginalXy d«v«lopcd bccauitt of dlii«tla£action ylth 
machanmtical programning tachnlquaai tha optimality critarion appiroach initially 
railed on intuitive optimality criteria auch aa conatant atraaa-ratio and 
uniform a train-energy denalty condltlona. More recently, optimality criteria 
(and aaaoclated re-deelgn equationa) have been derived from the Kuhn-Tucker 
condltlona (see, e.g. , Refi 9) for a conatrained minimization problem. 

The optimality criterion approach aeema especially well-suited to problems 
with a large number of design variables. Since our proposed design problems 
will have a moderate number of variables and since deriving efficient re-design 
equations for our problems is not Immediately straightforward, we initially 
prefer the mathematical programming approach over the optimality-criterion 
approach. 

A structural optimization computer program, called CONMIN, is available 
from NASA, It is this program that la used in our present work. ‘CONMIN 
is based on the mathematical nonlinear programming method of feasible dlrectiona. 

Nevertheless, if CON>flN proves to be too expensive computationally (or in 
any other way unsuccessful for our work) we may turn to other approaches. 
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2. BACKGROUND 

2,1 Formulation of Protvlem 

Becauae numarlcally-baaad optimization la bast carried out with diacrete 
variablea» tZ>e finite element technique atanda aa the moat logical choice for 
the blade model, A recent reaearch project (Kef. 10) haa reaulted in a finita^ 
element computer program that ia Ideally aulted to the work propoaed here. The 
program allowa for tapered » twlated finite elements in a rotating environment. 

The existing code can calculate natural frequencies, (with or without aero- 
dynamic terms) and forced response. 

Another Important aspect of the rotor blade optimization problem la the 
selection of the optimality criteria and constraints to be Imposed. Our 
design problem has certain features which are unusual compared to typical 
problems occurring In the structural optimization literature. There are 
basically three categories of criteria. In the first class, one would 
minimize weight given constraints on the natural frequencies (i.e. frequency 
"windows*')’ In this case, a constraint on rotary inertia is also implied since 
a rotor must have sufficient Inertia to autorotate. The advantage of this 
approach is that It is directly related to the physical realities of design, 

The disadvantage, however, ia that the first guess will probably not be feasible 
(that is will not have frequencies that fall in the "windows"), This can be a 
stumbling block to convergence. A second type of criteria is one in which the 
objective is to minimize the discrepancies between desired frequencies and 
actual frequencies. The constraint then becomes a window on autorotatlonal 
inertia. Although this avoids unfeasible solutions, it does not directly 
minimize weight (although weight is limited by the autorotatlonal constraint) . 

An objective function could be constructed that combined blass mass and frequency 
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plactiiMint, but th« relatlv* weightings of the two component? it not obvious. 

The third category of constraint is to minimise vibrations directly without 
r«tgard to frequency placement. Although this appears on the surface to be 
the perfect solution* there are problems. First* calculation of vibrations 
is an order-of-magnitude more difficult than the calculation of frequencies. 
Second* past efforts at this have resulted in strange designs* incompatible 
with standard helicopter practice. Third* there is still the problem of the 
weight-vibration trade-off. In this work* we intend to concentrate on the 
first two categories with some attention to the third. 

Another type of constraint Involved in the problem is the limitation on 
structural properties. The blade planform* airfoil, and twist art chosen 
by the aerodynamicist on the basis of performance. The structural engineer 
must choose his design to fit in the aerodynamic envelope given.. There are 
five structural parameters to be chosen: 1) flapping stiffness* 2) inplane 
stiffness, 3) torsional stiffness, 4) mass, and 5) torsional moment of 
inertia. In practice, these cannot be chosen completely Independently. 

Figure 1 shows the envelope of a jtyplcal blade section. All stiffness is 
assumed to reside in a box-beam of dimension b x h with thicknesses s and t. 

This beam is placed as far forward as possible (to keep the elastic axis 

near the 1/4 chord). Mass properties are due to the box-beam, skin* honeycomb* 
and two lumped masses. The lumped mass in the tip is typical of rotor blades 

and is used to keep the mass center forward of the aerodynamic center. A 

second lumped mass is included to allow independent choice of mass and mass- 
moment. The constraints of this construction are clear and are listed on the 
figure. 

In addition, there are minimum constraints on b* h, s, t to hold centifugal 
loads and to remain within manufacturable limits. For example a simple minimum 
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Figure 1. Typical Blade Cross-Section 
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contttiiint on ar«» could coma ftom tha cantrifugal eonac^aint (not conaidarlng 
banding atraia). Thtia, If ia tha naxinua atraaa and if f ia a aafaty factor 
than 

<0.) > £ I S (M.)!J*t,yA, 

■ i-J+1 ^ ^ ^ 


A, > 


n 

£( 2 

i"1+l 






Of couraa* whan wa antav tha vibratory-caaponaa phaaa of tha work (diacuaaad 
later) banding atraaaaa will ba included. 

Similarly, our early work Unich naglacta flutter boundariaa) can never- 
thalasa include flutter criteria In a almpllflad manner. Flrat, wa can 
chooae frequency placement auch that no coaleacence occura between flap-lag, 
flap-toraion, or lag-toralon. Second, we can constrain the five parametara 
in Figure 1 such that the mass center is always forward of tha 1/4-chord, a 
common design practice to prevent torsion-flutter In rotor blades. 

For specific examples, a smaller space of design variables may be used. 

For example, in the general case one could force the elastic axis to the 
quarter chord by choosing b «■ C/2 - 2h; or one could choose to allow no lumped 
mass in the blade interior (d ■ 0). This would reduce the number of design 
variables fi^om 5 to 4 or 3. If one considers flapping deflection only, then 
all variables except t and d may be fixed with no loss of generality. (This 
Is equivalent to simply using El and mass as variables.) Similarly, for 
Inplanc one can consider only s and d; and for torsion, one can consider only 
b and a. Thus, for each of the three uncoupled cases one has three possibilities: 
1) vary stiffness only, 2) vary inertia only, 3) vary stiffness and inertia. 
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2.2 Dtfinltlon of Tasks 

Ths worki as propostd hara, consists of four prliiiai^ phasas; 

1) solution of tha invarsa problan - givan a sat of fraquancias find tha aass 
and stiffnass distribution; 2) solution for optinun forcad rasponsa - givan 
a loading condition# mininiza blada strassas and hub raactions# sub,iact to 
waight rastrictions; 3) inprovanonts to includa flap-lag*torsion coupling; 
and 4) application to existing rotor blades - givan an existing blade# 
improve its dyitamic characteristics by realistic structural modification * 

These phasas are described in detail below. In each phase# work can be 
performed with ncndiroansional aquations* Mondimans ional quantities to be 
considered in optimization are p (mass distribution)# (flapping stiffness)# 
n_ (inplane stiffness), (torsional stiffness)# and e (c,g, - e,a. offset). 

m O 

The nondimensional parameters to be included as relatively fixed are 
Y (lock number)# 0 (twist), K (stiffness of control system), 6 (radius of 
^^yration about c.g.) and ^ (aspect ratio). 

Aerodynamics assumptions and other constraints are described in the 
individual tasks indicated below, 

2.3 Finite -Element Model 

Although tapered# twisted elements are within our capabilities# we 

introduce here a simpler case which is also of value, The stiffnesses Gd# 

El EX are assumed ta be constant along the length of the element. The 
zz yy G. 

mass of the element is assumed to be evenly distributed on the two nodes. 

Let the deflection of an element in the y and z directions at a distance 
X be denoted as w(x)’ and 'v(x) # for which the displacement models are assumed 
to be polynomials of third degree. The expressions are given as 


A 


Dimensional cases will also be considered# however 


# 


to keep insight. 
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w()i) - '4(2*’ - 31*^ + 1^) + - 2x*> 

i V 

- 4 A^x) - ^(x^ - ixh 


w(x) - ^(2x^ « 3Ax^ + A^) + ^(3Ax^ - 2x^) 
A-^ A^ 

- ^(x^ - 2Ax^ + A^x) - ^(x^ - 

r r 


whera v^, and Vg teprasant the banding dagraaa of fraadom in tha 
zx plana and u^i M^p and rapraaant tha banding dagraaa of fraadom 
in tha yx plana 


( 1 ) 


Tha strain anargy dua to banding daformation can ba axpraasad as 


A El .2 2 .2 .2 EI„ .2 2 

^0 t-2“<^) + EX ^ ^ Idx 

° ix y* 3» 3x^ ^ 3* 


(li) 


( 111 ) 


The potential anargy in tanaion from tha centrifugal force fialdf 
which is equivalent to the negative of kinetic energy due to radial 
displacement, is given by 


- T. 


0 T 


3wv2 /Svv2, . 

+ (gj) JdX 


where T, tension force, is assumed to be constant along each element. 
The kinetic energy due to inplane displacement is given by 


T^ - 1/2 / mvVdx 

N O 


which is equivalent to a u « -T 


-u 
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Meanwhile^ the pretwlet angle (|)(x)» and the toriional deformation 6(x) 
are aaaumed to be polynomlala of firat degree, and can be expreaaed aa 

<J><x) - -I) + 


e(x) - VjU - 1) + 

where ({i^j (|>2 repreaent the pretwiat angle at node 1 and 2, and v^, ''lO 
repreaent the elaatlc toralonaX degree of freedom at each end. 

(iv) The toraional energy » due to elaatlc deformatlona and centrifugal 
terms, can be expressed as 

u - kaj I <<||' + 0')^ + I ka^ j + 6')^ + | Gj6'^}dx 
where ka^ - ffz^dydz - 1^^ 
ka 2 » ffy^dydz - 

(v) The '*torslon-rotation” energy under the effect of rotation is given by 
t 2 

T^ ■ - (km^ - km^Htf) + 0)^dx 

2 2 

where km^, km 2 are mass moment of inertia which can be expressed as 
km^ *• // pr^dydz - 
kra? - // py^dydz “ pi 

^ Z 2H 
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Total diplacamant anargy now cPn ba uaad to fotw tha atlffnaaa matrix from 

— 1 T 

u ■ Y 

whara u ia the vector of nodal dlaplacamanta, In tha ordar aa u^, u^, 

Ug, U2» Mjf u^» Ug, My [k] la tha aXamantaX atlffnaaa matrix of ordar 10. 
(vi) The maaa matrix will ba obtained by the kinetic anargy of an element i 
which is given by 


/(B (iv.2 m .Jw.2 ^ / B^v a 2 


2 'at' 


2 *22 'atax^ 




Written in matrix form^ the kinetic energy can be expressed as 
f - I J [m)u 

where [m] is the mass matrix. 


13- 


ORrQINAl PAQg m 
OF POOR QUALITY. 


3. STATEMEHT OF WORK 


3,1 Phase 1 - The Inverse Problem 

The first phase of the research effort will Involve solution of an Inverse 
problem - that is, find the system structure given the eigenvalues. The 
solution to the inverse problem is difficult from a mathematical point of view, 
The Inverse problem that will be treated in this phase will be done in 
successive steps. First, a nonrotating cantilever beam will be used to study 
bending frequencies. Flapping will be considered ill the first beginning* 

Then, in-plane and torsion will be added. No twist or offset will be Included 
so that the two bending frequencies and the torsion frequency will be uncoupled 
A conventional optimal design formulation of one of the problems might be 
<1) minimize total weight 
n 

w - 2 jJ.Ar. 
i-1 ^ ^ 

With respect to design variables p and ri, subject to constraints 
E r^y^Ar^ -I 

- e < w. < + e 

i 11 


P 

n 


i ^ %in 


min 


> h > n 


max 


and w^ 


T 

X^KX 


T 

x^mx 


i 

i 


where K and m are structural stiffness and mass matrices with the rigid 
degrees of freedom removed, x is the eigenvector corresponding to the natural 
frequency, and w^ is the prescribed '.latural frequency. 
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The solution of the above optimization procedure is effected by the 
optimization program^ CONMIN. The efficiency of that program is greatly 
improved when the derivative of the constraints (with respect to the parameters) 
is known. Thus, we roust consider the derivative of frequency, with respect to 
the parameters : 




T 

X. in Xj 

J J 


3.2 Phase 2 > Forced Response 

In phase two, we will solve for blade structures that minimize forced 
response. In this step a very simplified model will be used in that we will 
simply ap| 5 ly a given forcing function. We will not take into account that 
the blade motions themselves may affect this loading, except that aerodynamic 
damping will be Included. We will assume a very simple lift distribution, such 
as u - Vl - r^ , that oscillates with integer-multiple frequencies 

Lift « E w u(r)e^“'^ 

— i ^ 


We will then foirmulate the problem as: 

A) Minimize (e {root shear] 4* |root moment]); that is, minimize 


n 

Z (w.u. 
l-X ^ ^ 




Ar^d + er^) 


given constraints a 


max 


I y ^ r Ar 
i-1 * 


“critical’ 

- 1. Pi > y 


min 


The forcing functions will be weighted probably as w^ ■ (.3)^ or some similar 
decreasing weightings. The a will include tension stress plus the oscillatory 
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stress* One of the major purposes of this will be to decide if minimizing huh 
loads iK equivalent to spacing frequencies half-way between integers* If this 
is so» then future optimizations would get by without a complicated loads 
calculation* There is some reason to suspect that this is so* Conceptually » 
the hub loads could be expressed as 


2 2 

i»j 1-w^/j^ 


2 2 

Thus, minimizing loads is similar to maximizing (1-w ) (1-w /4)*** 

2 2 2 2 

(1-w /n )* The simple factor (1-w )(l-w /4) has a maximum at w ■ 1*58 and 
2 2 

the factor (1-w /4)(l-w /9) has a maximum at w - 2.56. Thus, at (or slightly 
above) these points may well be optimum. 

3.3 Phase 3 - Improvements to Model 

In phase three, two parallel efforts are planned. In the first one, the 
blade structural model will be expanded so as to iticlude twist and c.g. offset. 
These will couple flap, inplane, and torsion. Therefore, uncoupled analyses 
will no longer be adequate. We already have computer programs that do this 
(Ref. 10) but the coupling will mean a more difficult convergence task for 
the optimization programs. In this phase, c.g. offset will be Included as a 
variable parameter. This implies that *'no flutter" must be an added 
constraint (Ref. 13). 

A parallel effort in Phase 3 will be the Improvement of the aerodynamic 
model to Include aerodynamic coupling between flap, lag, and torsion. The 
extent of this work will be partially determined by the results of Phase 2. 

Phase 4 - Application to Existing Designs 

Finally, these two parallel improvements (structural coupling and aerodynamic 
coupling) will be combined in a general program. Phase 4 will be to apply this 
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prograa to axlating rotor-blads daalgni to datannlna what atructural 
roodlHcations would ba taconmandad • Tha £act that our programa alraady agraa 
wall with axlatlng raaulta conatitutaa a atrong atarting point for atructural 
improvemanta* 
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4. ILLUSTRATIVE EXAMPLES OF STRUCTURAL OPTIMIZATION 

4.1 Cantilever Beana with Given Frequency 

Some simple examples will be examined and discussed before the utilization 
of the program CONMXN. In each case the results will be compared with those 
obtained by previous researchers, If It Is aviillable. 

The first limiting example Is the problem of determining the optimal design 
of an elastic cantilever beam, such that with a specific natural frequency, the 
weight of the structure attains the minimum value. 

We start with a uniform beam, modeled by ten elements, with a given 

2 3 

length of 10 Inches, E 1.0 Ib-ln , El * 10 lb-ln» density * 0*042 Ib/ln , 

and a specified first lowest natural frequency ■ 0.6489 rad/sec. We obtain 

the final stiffness profile shown In Figure 2, Figure 2 Is the present 

result with ten elements. 

4.2 Cantilever Beam with Given Weight 

A related problem has also been treated by Niels Olhoff [11]. He seeks 
the design of a cantilever beam that yields a maximum value of a particular 
higher natural frequency w^ (i.c., of specified order, n) with the volume 
and length of the beam specified. His work Is the dual problem of the example 
shown In Figure 2. Optimization with respect to the fundamental frequency 
under the constraint of volume Is similar to the one of minimizing weight 
(or volume) under the constraint of specified natural frequency. 

Figure 3 gives the profile of the optimal cantilever for n ■ 1 by Olhoff. 
One can see that the shapes In Figures 2 and 3 are very similar. 

4.3 Cantilever with Tip Mass 

Another example problem Is to minimize the weight of a cantilever carrying 
a mass at the tip, subject to the constraint that the fundamental natural 
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frcquency muit be than or aqual to a apacifiad valua. Tha problaat waa 

originally fomulatad by Turnar [12]. 

Kahn and Wlllmart [14] uaad an optimality critarlon nathod to aolva 

tumar'a problan* In this axamplop four finita alamanta ara uaad, vith tha 

areaa of each aa tha daaign variablaa, aa illuatrated in Flgura 4. Tha 

apacifiad natural fraquancy ia 17*752 rad/aac. Tha othar Initial data ara 

Modulua of alaaticity ■ 10*3 x 10^ pai 

-4 2 1 4 

Maaa danalty « 2.5 x 10 Ib-a /in 

Radlua of gyration (A^> • 2.0 

Radius of gyration (A 2 ) • 1*5 

Radius of gyration (A^) ■ 1 

Radius of gyration (A^) *0.5 

Concentrated roaart ■ 1 Ib-a /ln 

Length of each element ■ 60 in 

2 

where I ■ A (radius of gyration) . 
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Ths r«tults of th« optinlzAtion ar* ahowii in Tabla 1. Tha faaalblo 
•tatting daaign ia daactibad by ■ 200, A 2 " 150, A^ ■ 60 and A^ ■ 35. 


Table 1 


Compariaon of Cantilever Beam with Concentrated Maaa 



Ref. [12] 

Ref. [13] 

This Paper 

Iteration 


i3 

10 


136,81 

136.63 

134.60 

^2 

118.73 

118.7 

116.58 

^3 

83.591 

83.586 

82.734 

*4 

34.427 

34.608 

34.898 

Weight 

2243.0 

2242.9 

2214.41 


It can be aeen that excellent results have been obtained using the present 


CONMXM optimization program. 
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5. PRELIHXWAM CALCULATION FOR ROTORS 

5.1 Wind Tucbln* Blade 

The next example ie the optimization of a wind turbine rotor blade at 
30 rpm* A ten-element model ie ueed. Only the flapping ie conaldered. The 
area moment of inertia » end the lumped weight of each element ere taken 

o 

as the design variables. Young's Modulus » E « 0.2 x 10 Ib-ln, and 

3 

density *• 0.0334 Ib/in are assumed to be constant. Blade radius » R " 750 inches. 

Table 2 shows the profile of moment of Inertia and the distribution of 
added weight iot the initial and final configurations. 

The final profile of the area moment of inertia along the blade is 
similar as the one in the previous example. The lumped mass is concentrated 
at the tip of the blade as might be expected given a minimization of weight 
with a fixed moment of inertia. That is, the moment of inertia remains at 
the minimum value, as expected. We also note that most of the lumt'^ed mass 
is removed so that only mass necessary to the stiffness elements (or necessary 
for the autorotation constraint) is maintained. 

5.2 Other Considerations 

An important aspect of the optimization problem is the existence (or lack 
of it ) of a feasible solution. A "feasible solution" is defined as any set of 
design variables that satisfy the constraints (whether or not that particular 
solution is an optimum). It is possible that, if the problem is poorly 
formulated, that no feasible solution exists. What is more often the case, 
however, is that there are feasible solutions f ^t that the optimization scheme 
may not be able to find them. Thus, it is advantageous to have a feasible 
initial guess so that one is assured that at least a local optimum is possible. 
For example, Table 2 illustrates that the first guess is feasible 
(m^ > 2.82 per/rev) . Here we found that CONMIN was able to move from this 
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first gusss through th« tpsct ffsiiblt iolutions* In othsr caitt* hov«v«r, 
whtn the first gutss is not tsssibls vs hsvs found that CONHIN is not sbls to 
rssch A solution* In such cssss, ons must sdd or rsmovs soiss vsight (or add 
or rsstovt soma EX) from ths first gusss to movs into ths fsasibls spacs. From 
thsrs, optimisation is obtainsd* 

For sxampisy Tabls 3 prsssnts data for ths asms wind turbins as in Tabls 2, 
but ths constraint on ths first natural frsqusncy has bssn lowsrsd to rsmovs 
It from ths dangerous 3/rsv range • This implies that the first gusss in Tabls 2 
is no longer feasible* In order to overcome thiSi a lumped mass is added to 
station 9 (225*4 vs 49.50)* This lowers w^ below 2.621 rev but also lowers 
W 2 to 0.25/rev. This could be alleviated in one of two ways: 1) move the 
mass to the node of the second mode, or 2) simply widen the W 2 window. We have 
done the latter. It is interesting that the added weight Is ultimately 
rearrangea to other places and other weight removed such that the new design 
Is no heavisi: than the optimum In Tabl. 2. Furthermora, Wj 1. ral.ad to 8.57 
SO that the ’’widened window” had no effect on the solution. 

5.3 Uelicoptet Blade 

The design and analysis of a helicopter blade is discussed. Similar to 

Section 3, only flapping is considered and a ten-e^i^i^ent model is used. Density 

—3 3 

is constant along the blade and equal to 0.18 x 10 slugs/ in . Young’s 

7 2 2 

Modulus is equal fco 0.49 x 10 lb-in at the root equal to 0.585 x 10 lb-in 

elsewhere. Blade radius is equal to 193 Inches. Results are given In Table 

4 and 5. 

In Table 4, w^ is in the desired range but W 2 Is too small. Furthermore, 
the autorotatlonal inertia is larger than necessary. In this case, the CONMIN 
program is able to remove mass and stiffness in such a way to raise W 2 and 
lower w^. The mlnimuni I set as a constraint (0.4) is reached at e ary point 


Table 3. Wind Turbine vlth Added Mass 
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except the root. The root remelned high to keep ^ 1.05. The new blede le 
one- third the original mase. In Table 5, a atlffer initial gueea ia uaed and 
the frequency la restrained from significant decrease. Furthermore, 
near its maximum value. In this case, the program CONMIN would like to 
decrease El and m, but any removal of material could lower beyond its lower 
bound of 1.2A/rev . To counter this, the optimization scheme adds El near the 
root (to maintain > 1.24/rev). Furthermore, the lumped mass necessary 
to maintain ai^torotational constraint is moved slightly inboard to have less 
effect on w. (keep it high) but more effect on w» (keep it low). 
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6. SUMMAJIY 


Thus far we have completed the following objectlvea: 

1) Define the optimization problem and obtain physically reasonable 
constraints. 

2) Make CONMIN operational on our own computers and successfully apply 
It to beam problems with known solutions. 

3) Apply optimization procedures to typical rotor configurations 
for flapping. 

What remains In the first year Is 

1) Try many more flapping optimization runs to determine the effects 

of mass versus EX vard.atlons, to study the limits on range of frequency 
placement, to study the limit on number of frequencies placed, to 
determine to what extent a feasible first guess Is necessary. 

2) Determine the feasibility of using frequency-placement as an objective 
function rather than as a constraint 

3) Try Inplane and torsion optimization. 

After this, we Intend to follow the statement of work as planned. 
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8. NOTATION 

pXopa of lilt curve, rad 

size of tip maaa, m 

width of box beam, m 

blade chord, m 

width of lumped mass, m 

weighting function 

safety factor 

height of box beam, m 

nondlmenslonal mass/ unit length, m R/m 

Gd/m^R^ 

EI^/m^R^ 

lock number, pacR^/1 

control system stiffness, N-m/rad 

K/m^R^ 

blade twist, rad 

nondlmenslonal c.g. elastic-axis offset (offset/R) 

radius of gyration 'Vl^/l-R^ 

inverse aspect ratio c/r 

density of air 

small parameter 

Qt 

rotating speed 
maximum stress 
mass of skin kg/m 

3 

density of lumped mass, box beam, honeycomb kg/ra 


I 
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2 

nic«teary Inareia for «ucorotatlon» kg-m 
toralotial maia momant of Inartla (about c.g.) 
par unit langth* kg-m 

2 

flapping atiffnaaai N-m 

2 

Inplana atlffnaas, N-m 

2 

torilonal atlffnaaBi N-m 

maaB/unit length | kg/m 

2 

reference mass, I/R , kg 
x/R 

length of element, m 
blade radlua, m 
thicknesses of box beam, m 
lift distribution, N/m 
displacements, m 
weighting functions 
length along blade, m 


